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We observe a narrow enhancement near 2m p in the in- 
variant mass spectrum of pp pairs from radiative J/ip — + 7pp 
decays. No similar structure is seen in J/ip — » Tx°pp decays. 
The results are based on an analysis of a 58 million event 
sample of J/ip decays accumulated with the BESII detector 
at the Beijing electron-positron collider. The enhancement 
can be fit with either an S- or P-wave Breit Wigner reso- 



nance function. In the case of the S-wave fit, the peak mass 
is below 2m p at M = 1859_+q ( stat )-~25 ( s y s ) MeV/c 2 and the 
total width is F < 30 MeV/c 2 at the 90 percent confidence 
level. These mass and width values are not consistent with 
the properties of any known particle. 
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There is an accumulation of evidence for anomalous 
behavior in the proton-antiproton (pp) system very near 
the M P p = 2m p mass threshold. The observed cross sec- 
tion for e + e~ — > hadrons has a narrow dip-like struc- 
ture at a center of mass energy of y/s ~ 2m p c 2 [1]. The 
proton's time- like magnetic form-factor, determined from 
high statistics measurements of the pp — > e + e~ annihi- 
lation process, exhibits a very steep fall-off just above 
the pp mass threshold [2]. These data arc suggestive of 
a narrow, S- wave triplet pp resonance with J PC = 1~~ 
and mass near 2m p . In studies of p annihilations at rest 
in deuterium, anomalies in the charged pion momentum 
spectrum from pd — > Tr~n°p and 7r + 7r~n reactions [3] 
and the proton spectrum from pd — > p27r + 37r~ [4] have 
been interpreted as effects of narrow, below-threshold res- 
onances. There are no well established mesons that could 
be associated with such states. The proximity in mass 
to 2m p is suggestive of nucleon-antinuclcon (NN) bound 
states, an idea that has a long history. In 1949, Fermi 
and Yang [5] proposed that the pion was a tightly bound 
NN state. Nambu and Jona-Lasinio [6] expanded on this 
in 1961 with a model based on chiral symmetry that has, 
in addition to a low-mass pion, a scalar pp composite 
state with mass equal to 2m p . Although these ideas have 
been superseded by the quark model [7], the possibility 
of bound NN states with mass near 2m p , generally re- 
ferred to as baryonium, continues to be considered [8]. 
Recently Belle has reported observations of the decays 
B+ -> K+pp [9] and B° -> D°pp [10]. In both pro- 
cesses there are enhancements in the pp invariant mass 
distributions near M pp ~ 2m p . An investigation of low 
mass pp systems with different quantum numbers may 
help clarify the situation. 

In this letter we report a study of the low mass pp 
pairs produced via radiative decays in a sample of 58 
million J/ip events accumulated in the upgraded Beijing 
Spectrometer (BESII) located at the Beijing Electron- 
Positron Collider (BEPC) at the Beijing Institute of High 
Energy Physics. For this reaction, charge-parity conser- 
vation insures that the pp system has C = +1. 

BESII is a large solid-angle magnetic spectrometer 
that is described in detail in rcf. [11]. Charged parti- 
cle momenta are determined with a resolution of u p /p — 

1.7%\Jl+p 2 (GeV 2 ) in a 40-1 ayer cylindrical drift cham- 
ber. Particle identification is accomplished by specific 
ionization (dE/dx) measurements in the drift chamber 
and time-of-flight (TOF) measurements in a barrel-like 
array of 48 scintillation counters. The dE/dx resolution 
is (JdE/dx = 8.4%; the TOF resolution is otof = 180 ps; 
both systems independently provide more than 3cr sepa- 
ration of protons from any other charged particle species 
for the entire momentum range considered in this exper- 
iment. Radially outside of the time-of-flight counters is 
a 12-radiation-length barrel shower counter (BSC) com- 
prised of gas proportional tubes interleaved with lead 



sheets. The BSC measures the energies and directions 
of photons with resolutions of a E /E ~ 22%/^/E(GeV), 
(Jcj, — 4.5 mrad, and ag = 7.9 mrad. The iron flux return 
of the magnet is instrumented with three double layers 
of counters that are used to identify muons. 

For this analysis we use events with a high energy 
gamma ray and two oppositely charged tracks each of 
which is well fitted to a helix originating near the inter- 
action point. Candidate 7's are associated with energy 
clusters in the BSC that have less than 80% of their total 
energy in any one readout layer and do not match the ex- 
trapolated position of any charged track. Since antipro- 
tons that stop in the material of the TOF or BSC can pro- 
duce annihilation products that are reconstructed else- 
where in the detector as 7 rays, no restrictions are placed 
on the total number of neutral clusters in the event. We 
use charged tracks and 7's that are within the polar angle 
region | cos0| < 0.8 and reject events where both tracks 
arc identified as muons, or produce high energy show- 
ers in the BSC that are characteristic of electrons. The 
dE/dx information is used to form particle identification 
confidence levels V pid , where i denotes ir, K and p. We 
require that both charged tracks have V pid > V pid and 
V pid > Vp id - A study based on a kinematically selected 
sample of J/ip — > K* ± K T —> K + K~ir° events indicate 
that the probability for a charged kaon to satisfy this 
requirement is less than 1% per track. 

We subject the surviving events to four-constraint 
kinematic fits to the hypotheses J/ip — > jpp and J/ip — > 
■yK + K~. For events with more that one 7, we select 
the 7 that has the highest fit confidence level. We select 
events that have fit confidence level CL lpp > 0.05 and 
reject events that have CL yK + K - > CL lpp . 
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FIG. 1. The pp invariant mass distribution for the 
J/ip — > 7pp-enriched event sample. 



Figure 1 shows the pp invariant mass distribution for 
surviving events. The distribution has a peak near 
M pp — 2.98 GeV/c 2 that is consistent in mass, width, 
and yield with expectations for J/ip — > jr) c , r\ c — > pp [12], 
a broad enhancement around M P p ~ 2.2 GeV/c 2 , and a 
narrow, low-mass peak at the pp mass threshold that is 
the subject of this Letter. 
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Backgrounds from processes involving charged parti- 
cles that are not protons and antiprotons are negligibly 
small. In addition to being well separated from other 
charged particles by the dE/dx measurements and the 
kinematic fit, the protons and antiprotons from the low 
M p p region stop in the TOF counters and, thus, have very 
characteristic BSC responses: protons do not produce 
any matching signals in the BSC while secondary parti- 
cles from antiproton annihilation usually produce large 
signals. This asymmetric behavior is quite distinct from 
that for K + K~, n + n~ or e+e~ pairs, where the posi- 
tive and negative tracks produce similar, non-zero BSC 
responses. The observed BSC energy distributions for 
the selected J/iji — ► "fpp events with M pp < 1.9 GeV/c 2 
closely match expectations for protons and antiprotons 
and show no evidence for contamination from other par- 
ticle species. 

There is, however, a large background from J/ip — ► 
ir°pp events with an asymmetric ir° — > 77 decay where 
one of the photons has most of the n°'s energy. This 
is studied using a sample of J/ip — > n°pp decays re- 
constructed from the same data sample. For these, we 
select events with oppositely charged tracks that are 
identified as protons and with two or more photons, 
apply a four-constraint kinematic fit to the hypothesis 
J/ip — ► 77PP, and require CL 11V p > 0.005. For events 
with more than two 7's, we select the 7 pair that pro- 
duces the best fit. In the M 77 distribution of the se- 
lected events there is a distinct ir° signal; we require 
|M 77 - M n o\ < 0.03 GeV/c 2 (±2<r). The distribution 
of events vs. M pp - — 2m p near the M pp = 2m p thresh- 
old, shown in Fig. 2(a), is reasonably well described by a 
function of the form f hhg (S) = N{5^ 2 + M 3 / 2 +a 2 <* 5/2 ), 
where 5 = M pp — 2m p and the shape parameters a\ and 
<i2 arc determined from a fit to simulated MC events that 
were generated uniformly in phase space. This is shown 
in the figure as a smooth curve. There is no indication of 
a narrow peak at low pp invariant masses. Monte Carlo 
simulations of other J/ ip decay processes with final-state 
pp pairs indicate that backgrounds from processes other 
than J/ip — > 7r°pp are negligibly small. 

The M pp — 2m p distribution for the Tr°pp phase-space 
MC events that pass the jpp selection is shown in 
Fig. 2(b). There is no clustering at threshold; the smooth 
curve is the result of a fit to /bkg(^) with the same shape 
parameter values. 

In BESII, the detection efficiency for protons and 
antiprotons falls sharply for three-momenta below 
0.4 GeV/c. This produces a mass dependence in the ex- 
perimental acceptance near M pp ~ 2m p for J/ip — > -fpp 
and ir°pp. For both processes, when M pp is very near 
2m pi the p and p both have three- momenta very near 
0.5 GeV/c and are well detected. For increasing pp 
masses, more asymmetric energy sharing is possible and 
the acceptance decreases until M pp ~ 2.0 GeV/c 2 , where 
it is ~ 0.65 of its value at M™ = 2m v . 
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FIG. 2. The M P p — 2m p distribution for (a) selected 
J/ip — » 7r°pp decays and (b) MC J/ip — > TT°pp events that 
satisfy the 7pp selection criteria. The smooth curves are the 
result of fits described in the text. 

Figure 3(a) shows the threshold region for the selected 
J/ip — ► jpp events. The dotted curve in the figure indi- 
cates how the acceptance varies with invariant mass. The 
solid curve shows the result of a fit using an acceptance- 
weighted S'-wave Breit-Wigner (BW) function [13] to rep- 
resent the low-mass enhancement plus fbk s (S) to repre- 
sent the background. The mass and width of the BW 
signal function are allowed to vary and the shape pa- 
rameters of fbkg(S) are fixed at the values derived from 
the fit to the TT°pp phase-space MC sample [14]. This 
fit yields 928 ± 57 events in the BW function with a 
peak mass of M = 1859 +^ MeV/c 2 and a full width 
of T = 0± 21 MeV/c 2 [15]. Here the errors are statistical 
only: those for the event yield and the width are derived 
from the fit; the determination of the statistical errors 
for the mass is discussed below. The fit confidence level 
is 46.2% ( X 2 /d.o.f. = 56.3/56). 

Monte Carlo studies indicate that in the presence of 
background, the determination of the peak mass for a 
below-threshold resonance is more unreliable the further 
the peak position is below threshold. This produces an 
asymmetric distribution of mass input values that can 
produce our measured result. Moreover, the rms spread 
of these values increases for lower input masses, indi- 
cating that the statistical error returned by our mass 
fit underestimates the negative error. Because of this, 
we quote statistical errors for the mass that are derived 
from the rms spreads of fit results for an ensemble of MC 
experiments with different input mass values. 
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FIG. 3. (a) The near threshold M P p — 2m p distribution for 

the 7pp event sample. The dashed curve is the background 

function described in the text. The dotted curve indicates how 

the acceptance varies with pp invariant mass; the dashed curve 

shows the fitted background function, (b) The M pp - — 2m p 

distribution with events weighted by qa/q. 

Further evidence that the peak mass is below the 2m p 
threshold is provided in Fig. 3(b), which shows the M P p — 
2m p distribution when the kinematic threshold behavior 
is removed by weighting each event by qo/q, where q is 
the proton momentum in the pp restframe and qo is the 
value for M pp - = 2 GeV/c 2 . The sharp and monotonic 
increase at threshold that is observed in this weighted 
histogram can only occur for an S'-wave BW function 
when the peak mass is below 2m p . 

An S-wave pp system with even C-parity would cor- 
respond to a h pseudoscalar state. We also tried 
to fit the signal with a P-wave BW function, which 
would correspond to a ++ (pP) scalar state that oc- 
curs in some models [6,8]. This fit yields a peak mass 
M = 1876.4 ± 0.9 MeVc 2 , which is very nearly equal to 
2m p , and a very narrow total width: T = 4. 6 ± 1.8 MeVc 2 
(statistical errors only). The fit quality, x 2 /d.o.f. = 
59.0/56, is worse than that for the S-wave BW but still 
acceptable. A fit with a P-wave BW fails badly with 
X 2 /d.o.f. = 1405/56. 

In addition we tried fits that use known particle res- 
onances to represent the low-mass peak. There are 
two spin-zero resonances listed in the PDG tables in 
this mass region [16]: the ?7(1760) with A^(i76o) = 
1760 ± 11 MeV/c 2 and r^ (1760) = 60 ± 16 MeV, and 
the tt(1800) with M <1800) = 1801 ± 13 MeV/c 2 and 
1^(1800) = 210 ± 15 MeV. A fit with / bkg and an 
acceptance- weighted S'-wave BW function with mass and 
width fixed at the PDG values for the 77(1760) produces 
X 2 /d.o.f. = 323.4/58. A fit using a BW with the tt(1800) 
parameters is worse. 

For both the scalar or pseudoscalar case, the polar 



angle of the photon, 6L, would be distributed accord- 
ing to 1 + cos 2 6> 7 . Figure 4 shows the background- 
subtracted, acceptance-corrected | cos# 7 | distribution for 
events with M pp < 1.9 GeV and |cos6> 7 | < 0.8. Here 
we have subtracted the | cos 9^o | distribution from the 
ir°pp data sample, normalized to the area of /bkg(^) for 
M pp < 1.9 GeV/c 2 to account for background. The 
solid curve shows the result of a fit for 1 + cos 2 8 y to 
the I cos^l < 0.8 region; the dashed line shows the re- 
sult of a similar fit to sin 2 # 7 . Although the data are not 
precise enough to establish a l + cos 2 9 1 behavior, the dis- 
tribution is consistent with expectations for a radiative 
transition to a pseudoscalar or scalar meson [17]. 
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FIG. 4. The background-subtracted, acceptance-corrected 
I cos 6~f I distribution for J/ip — > 7pp-enriched events with 
M P p < 1.9 GeV/c 2 . The solid curve is a fit to a 1 + cos 2 # 7 
shape for the region |cos# 7 | < 0.8; the dashed curve is the 
result of a fit to sin 2 9-, . 



We evaluate systematic errors on the mass and width 
from changes observed in the fitted values for fits with dif- 
ferent bin sizes, with background shape parameters left 
as free parameters, different shapes for the acceptance 
variation, and different resolutions. The ensemble Monte 
Carlo studies mentioned above indicate that in the pres- 
ence of background, the determination of the parameters 
of a sub-threshold BW resonance can be biased. We in- 
clude the range of differences between input and output 
values seen in the MC study in the systematic errors. 

For the mass, we determine a systematic error of 
MeVc 2 . For the total width, we determine a 90% 
confidence level (CL) upper limit of T < 30 MeV/c 2 , 
where the limit includes the systematic error. 

Using a Monte-Carlo determined acceptance of 
23%, we determine a product of branching fractions 
B(J/ip -> 7V(1859))S(V(1859) -> pp) = (7.0 ± 
0.4(stat)^ ( ! ) 'g(syst)) x 10~ 5 , where the systematic error 
includes uncertainties in the acceptance (10%), the total 
number of J/ip decays in the data sample (5%), and the 
effects of changing the various inputs to the fit ( + _^2%)- 

In summary, we observe a strong, near-threshold en- 
hancement in the pp invariant mass distribution in the 
radiative decay process J/ip — ► jpp. No similar struc- 
ture is seen in J/ip — > Tr°pp decays. The structure has 
properties consistent with either a J PC = ^ or ++ 
quantum number assignment and cannot be attributed 
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to the effects of any known meson resonance. If inter- 
preted as a single h resonance, its peak mass is be- 
low the M P p = 2m p threshold at 1859j^g(stat) jt^syst) 
MeV/c 2 and its width is V < 30 MeV/c 2 at the 90% CL. 
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